Condensation transition of ultracold Bose gases with Rashba spin-orbit coupling 
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We stxidy the Bose-Einstein condensate phase transition of three dimensional ultracold bosons with 
isotropic Rashba spin-orbit coupling. Investigating the structure of Ginzburg-Landau free energy as a 
function of the condensate density, we show that the condensate phase transition is first order with a jump 
in the condensate density. We calculate the transition temperature and the jump in the condensate density 
at the transition for large spin-orbit coupling, where the transition temperature depends linearly on the 
density of particles. Finally, we discuss the feasibility of producing the phase transition experimentally. 
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The recent experimental realization of artificial non- 
Abelian gauge fields in ultracold atomic gases Q-lsl] has 
opened the prospect of exploring a wide variety of mod- 
els that are difficult to realize or investigate in other phys- 
ical systems. Of particular interest is simulated Rashba- 
Dresselhaus spin-orbit coupling for which the modi- 
fied single-particle spectrum leads to novel phases and phe- 
nomena. A notable feature of such systems is that the nor- 
mal phase is not kinematically forbidden at any nonzero 
temperature, unlike in usual three dimensional Bose gases 
where there is no normal state below the EEC transition 
temperature. On the other hand, interparticle interactions 
stabilize Bose condensates at low temperatures fll3]- The 
issue we investigate in this paper is the phase transition at 
a nonzero temperature between the normal and condensed 
states. 

The ground state of ultracold bosons with Rashba- 
Dresselhaus spin-orbit coupling is predicted, depending on 
the details of the interactions, to have two characteristic 
phases: a plane-wave state, which is a Bose-Einstein con- 
densate (BEC) in a single momentum state, or a striped 
state, which is a BEC of two opposite momenta |.8. -12ll . 
More exotic phases are predicted at nonzero tempera- 
tures (]3\ and in harmonic traps lfl4l - [l8!l . Atoms with 
isotropic in-plane Rashba spin-orbit coupling have circu- 
larly degenerate single particle energy minima; as a conse- 
quence non-interacting bosons with such coupling do not 
Bose condense in three dimensions, owing to the density 
of states being two-dimensional at low energy lEl. On 
the other hand, as mentioned, interparticle interactions al- 
low Bose condensation at nonzero temperatures. Thus the 
phase of the system, whether condensate or normal, is de- 
termined from energetics. 

We focus here on spatially homogeneous ultracold 
bosons in the three dimensions with isotropic in-plane 
Rashba spin-orbit coupling. Assuming a plane-wave con- 
densate at low temperatures, we explore the condensate 
phase transition by obtaining the Ginzburg-Landau free en- 
ergy as a function of the condensate density around the 
transition, within the Bogoliubov approximation includ- 
ing Hartree-Fock energies (henceforth referred to as BHF). 
Calculating how the coefficients of the Ginzburg-Landau 
free energy vary with temperature, we conclude that the 



transition between the condensed and normal phases is first 
order. For relatively large spin-orbit coupling, we estimate 
the transition temperature, which depends linearly on the 
density of the particles, n. With increase in the spin-orbit 
coupling, scattering length or particle density, it can be- 
come feasible to see the transition in realistic experimental 
setups. 

In usual Bose gases one finds a spurious first order 
phase transition within BHF, driven by order parameter 
fluctuations. Such fluctuations lead to a density of par- 
ticles excited out of the condensate near the transition, 

1/2 

~ n — auQ , which is non-analytic in the condensate frac- 
tion, Uq. The first order transition is removed when cor- 
rectly determining the critical behavior at the phase tran- 
sition [20]. The relevant momentum scale in this problem 
is ^JariQ where a is the s-wave scattering length. How- 
ever, Rashba spin-orbit coupling introduces a second scale, 
K, the strength of the spin-orbit coupling (see Eq. ([T]l be- 
low); we see as a consequence that the density of excited 
particles is analytic in Uq for ^ ang, where a is the 
interparticle s-wave scattering length. Whether a more ex- 
act treatment beyond mean field would make the transition 
second order remains a problem for the future. 

Hamiltonian. We consider a system of bosons in two 
hyperfine (pseudospin) states labelled a and b, with an 
isotropic in-plane Rashba-Dresselhaus spin-orbit coupling 
and an isotropic s-wave interaction, described by the 
Hamiltonian 
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As before ||6|], m is the atomic mass, V the volume of the 
system, k the spin-orbit coupling strength, taken to be pos- 
itive, and Op and bp are annihilation operators of particles 
of momentum p in the pseudospin states a and b. The 
(Ta; and ay are Pauli matrices between the internal states, 
and / is the two-by-two identity matrix. We assume an 
isotropic (constant) mean-field coupling g; extending the 
present theory to include the effects of renormalization of 
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the interaction rtllLl2lLl22ll is left as a future problem. Diag- 
onalization of the single particle terms in the Hamiltonian 
gives two single-particle dispersion branches 



e±{v) = {{p^±Kf +pl}/2m, 



(2) 



where p± = \Jpt + P^, with circularly degenerate ground 

states along {p±,Pz) = 0). A previous study 10] 
shows that the plane-wave state, in which the condensate 
is made of particles with a single momentum, is the pre- 
ferred ground state within the Bogoliubov approximation. 
We start from the plane -wave ground state with momen- 
tum K, = (k, 0, 0), and consider how the transition to the 
normal state takes place at finite temperatures. 

Normal state. We consider the free energies of the 
normal and condensed states as functions of temperature, 
T, and chemical potential, fi; at the phase transition, the 
chemical potentials and pressures of the normal and con- 
densate phases must be equal. In the normal state the quasi- 
particle dispersion relation within Hartree-Fock is 

^ , , (p± ± + pI 3 , , 

^±(P) = ^ ^^-/i + -5n(M), (3) 

where the total number of particles in the normal phase 
n(/i), which is a function of the chemical potential /x and 
the temperature T, is self-consistently determined from the 
number equation 

"(^) = ^E{/(^(p)) + /(^+(p))}' (4) 



where f(x) = l/{e^/'^ — 1). For a given total density n, 
the chemical potential approaches 3gn/2 as T — )• 0, and 
for nonzero temperatures, /i < 'ign/2. It is convenient to 
define the deviation of the chemical potential from its zero 
temperature value as 



A/i(^,r) 



(5) 



Since n(^) increases as a function of T, |A/i| also in- 
creases as T increases. 

Ginzburg-Landau free energy. We now determine, 
within BHF, the Ginzburg-Landau free energy as a function 
of the condensate density, no, around the condensate tran- 
sition,. Since the operator (a], — h\) /\/2 creates a particle 
in the plane -wave condensate, it is easier, as before |@|, to 
work in the following (— , +) basis: 



(6) 



The state created by V'-.k is macroscopically occupied. 
The effective Hamiltonian in the (— ,+) basis within 



BHF is 
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where no is the density of condensate particles, n_ and n+ 
are the densities of particles in (— ) and (-|-) states that are 
not in the condensate. The derivative of the free energy, 
T, no), with respect to no is 
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+ ^ E (V'^.pV'1.2.-p + ^-.p^-.2.-p) . (8) 
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As no — )• 0, one recovers the free energy of the normal 
phase n_,n+ — )• n{^)/2, and the last term approaches 
zero in this limit. In the following, we expand the right 
side of ([8]) as a function of ng to obtain the difference of the 
free energies T in the condensed and normal phases. The 



expansion is facilitated by using the single particle matrix 
Green's functions with anomalous components, 

G(q,ti -ts) = -i{T{^^{tr)^\{h))), (9) 
where the four-component spinor ^q(t) is 

(t), V'-,,^-q(i), V'+.K+qW, V'+,,c-q(i)) ■ 

(10) 

In terms of G, Eq. (jSj becomes 

dJ- f d?q ( 

+ G33(p,Z,) + ^(G2i(p,Z.) + Gi2(p,2.))j , (11) 
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where the v are the bosonic Matsubara frequencies. The 
matrix Green's function within BHF is tS, 
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(12) 



where 



^(q) = g^/2m - fi + g{2no + 2n_ + n+) 
S(q) = {2k + c\Y/2m - fi + g{no + n_ + 2n+) 
D(q) = i?(-q). (13) 

Since the condensate transition is characterized by the 
infrared structure of the Green's functions, in expanding 
in tlq we consider only the 1^ = component in ([TTT l. In 
evaluating Green's functions in (fTTll . we approximate the 
n_ and n+ in A, B, and L> through [\3\ by n(/i) /2, their 
value to leading order in no; including the Uq dependence 
of n± by solving the number equations self-consistently 
remains a task for the future. The right side of ([TTT l. which 
is analytic in griQ, can be expanded for small gno as 

-— = -/U + -gn{fi) + Xgno + Y{gnof + . . . , 
orio 2 



(14) 



where 



Y^,,^>.'-I^](^Y (,5, 

e,; = k'^ /2m, and Q!(x) and /3(a;) are dimensionless func- 
tions. Note that if one replaces g by Ana/m, the prefac- 
tor Am^g^T becomes 327r^a/A^, where a is the scattering 
length and A is the thermal wavelength. 

Since A/i < 0, we need only consider a{x) and /3(x) 
with negative arguments; by explicit calculation, one sees 
that a{x) and (3{x) are both positive there and monotoni- 
cally increasing functions of x, with the asymptotic forms 
as a; — )• 0~, 



a[x] ~ 
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0-19 or ^ 0.16 

> '• ^, (16) 



and approaching as x — )• — oo. The analytic form for a is 
derived from the 2Gii + G33 term in Eq. ([TTT l; the contri- 
bution of the G12 terms is small, and numerically changes 
the coefficient from 0.19 to 0.20. Figure [T] plots a{x) and 
fi{x), calculated numerically. 

Integrating (fT4] |. we obtain the Ginzburg-Landau free en- 
ergy up to third order in tiq: 



^(no) = -Fn - A/i no + X^nl + Y 
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-n, 



(17) 




FIG. 1: Functions a{x) (solid line) and /3(a;) (dashed line), 
Eq. l |15l >, for > x > —1. The inset shows the behavior of 
the functions in the lower right comer, where — a; > 1. 



where = J~{nQ = 0) is the free energy in the nor- 
mal phase. The coefficients of and n^ are both pos- 
itive, whereas for given /x > 0, the coefficient X of np 
is negative at low T, and decreases continuously with de- 
creasing temperature ['23^. This change in the coefficient 
of nQ drives a first order phase transition, since at suf- 
ficiently small T, the two conditions J- {no) = J-'^ and 
9J-"(no)/5no = become simultaneously satisfied; this 
occurs when X'^ = —16Y Afj,/3. At this temperature the 
system undergoes a transition to the condensed phase. Fig- 
ure |2] schematically shows how this transition takes place. 
The combination X"^ /(—Y Afi) = T monotonically de- 
creases with T, as long as X < 0. At the transition, Uq 
jumps from zero to |3X/4y(7| > on the condensate side. 




FIG. 2: The first order phase transition from the normal to con- 
densed phase atr = 16/3. . The lines show the Ginzburg-Landau 
free energy, measured with respect to the free energy of the nor- 
mal phase. The top line is for T > Tc where F < 16/3; the 
middle line is at Tc, where F = 16/3; and the bottom line is for 
T <Tc where F > 16/3. 



Transition temperature. Now we estimate the transition 
temperature assuming that the spin-orbit coupling strength 
is sufficiently large that 3> |A/i|. We will see that 
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this condition is obeyed by typical experimental parame- 
ters. Then, using the asymptotic forms [\M the condition 
for the transition, T = 16/3, becomes 



1 A/x 1 
2mKg Tj, 0.16 



1+0.20 X Imng-^ 
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(18) 



The left side depends only on the combination 
Afi{Tc)/{2mKgTc), which thus /S.^/ {2mKgTc) takes 
a constant universal value. Numerically solving ([Tsl l. 
choosing the solution which makes X < 0, we find 



-0.033. 



(19) 



2'mKgTc 

The shift in chemical potential can be written in 
terms of the normal state density n(/u) using the number 
equation dH. For |A^| <^ e^, and T/e^ ^ 1, the main 
contribution to dU comes from pj^ ~ k and ~ in 
^ /(C- (p)) /y ^ and we obtain 

n{^) « In i-Afi/T) . (20) 

zvr 

Then, using ( [T9l l and (|20] |. we obtain the transition temper- 
ature, for 2mKg ^ 1, in terms of the normal state density 
as 
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ln{2mKg)\ + C" 



(21) 



where C ~ 3.4. Thus the transition temperature depends 
linearly on the density, and as expected approaches zero 
as (jr — )■ 0, since there is no Bose condensate at nonzero 
temperature in the absence of interactions. 

The jump in the condensate density at the transition is 

^ . ^1^? , (22) 

n{fi) \ln{2mKg)\ + C 

approaching zero as (7 — t- 0, and increasing as g increases. 

Experimental feasibility. Finally we mention accessing 
the transition experimentally. For orientation we take typ- 
ical current experimental values from Ref. ijj], which re- 
alized a gas of ^''Rb with a mixture of Rashba and Dres- 
selhaus spin-orbit coupling, with coupling strength k ~ 
\/2vr/800nm and density of order n ~ 10^^/cm^. Ap- 
proximating the coupling g by iira/m, so that 2mK,g ~ 
0.69, we obtain T^, ~ 2nK with a jump in the conden- 
sate fraction at the transition, ng/n ~ 0.1. In addition, 
rje„ ~ 0.019 and A/i/e„ ~ 4.3 x lO"'' putting the sys- 
tem in a regime where the approximations leading to (|2T] | 
are valid. Future experiments could, depending on their 
specific configurations, be able to increase the transition 
temperature and the jump in the condensate density via in- 
creasing K, n, or a. 
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